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ABSTRACT 

Gamma-ray bursts (GRBs) are believed to originate from ultra-relativistic winds/fireballs to avoid the "com- 
pactness problem". However, the most energetic photons in GRBs may still suffer from 7 — 7 absorption leading 
to pair production in the winds/fireballs. We show here that in a wide range of model parameters, the resulting 
pairs may dominate those electrons associated with baryons. Later on, the pairs would be carried into a reverse 
shock so that a shocked pair-rich fireball may produce a strong flash at lower frequencies, i.e. in the IR band, in 
contrast with optical/UV emission from a pair-poor fireball. The IR emission would show a 5/2 spectral index due 
to strong self-absorption. Rapid responses to GRB triggers in the IR band would detect such strong flashes. The 
future detections of many IR flashes will infer that the rarity of prompt optical/UV emissions is in fact due to dust 
obscuration in the star formation regions. 

Subject headings: gamma-rays; bursts — radiation mechanisms: nonthermal — relativity 
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1. INTRODUCTION 

There is little doubt that GRBs arise from relativisti- 
cally expanding winds/fireballs, which also drive relativistic 
blast waves that responsible for later time emission at X- 
ray, optical and radio bands (see reviews of Piran 1999; van 
Paradijs, KouveUotou & Wijers, 2000; Cheng & Lu 2001; and 
Meszaros 2002). The observed GRB spectra are hard, with a 
significant fraction of the energy above the pair forma- 
tion energy threshold, so it is believable that there are pro- 
cesses of substantial pair production related to GRBs. One case 
of pair production in GRBs is that the GRB photons are 
back-scattered by the ambient medium and interact with other 
outgoing high-energy GRB photons, leading to pair pro- 
duction and then deposition of momentum into the external 
medium. Such a case has been considered by Madau & Thomp- 
son (2000); Thompson & Madau (2000); Dermer & Bottcher 
(2000); Madau, Blandford & Rees (2000); Meszaros , Ramirez- 
Ruiz & Rees (2001); Beloborodov (2002); Ramirez-Ruiz, Mac- 
Fadyen & Lazzati (2002), which could be called an "external 
case". In this paper, we consider the so-called "internal case", 
in which pairs are produced within the relativistic fireball due to 
7 — 7 absorption between high energy photons, and the formed 
pairs remain mixing in the fireball matter The pairs in the 
external medium (external case) will modify the usual prop- 
erties of the forward shock and then the later afterglow, while 
the pairs in the wind/fireball (internal case) affect the reverse 
shock that is responsible for prompt emission in soft bands dur- 
ing the bursting phase. Pilla & Loeb (1998) have also studied 
spectral signatures of GRBs by considering pair production in 
fireballs. It is possible that the pairs are so abundantly gen- 
erated that the fireball becomes optically thick, and the GRB 
photons undergo Compton multi-scattering before leaking out. 
Meszaros et al. (2002) discussed this case to account for the 
currently known X-ray flashes. 

A large Lorentz factor is required to avoid the "compactness 
problem" in GRBs (see, e.g., Piran 1999). However, the high- 
energy photons may still suffer from absorption, and the gener- 
ated pairs may dominate the electrons associated with baryons, 
and then change the conditions in the reverse shock. The pres- 



ence of a large number of pairs decreases the shared energy 
per lepton in the reverse shock, so that the reverse flash is soft- 
ened to lower frequencies. The reverse-shock model (Sari & Pi- 
ran 1999; Meszaros & Rees 1999) has been successful to inter- 
pret the strong optical flash detected in GRB 990123 (Akerlof 
et al. 1999). Here we show that a pair-rich reverse shock gives 
rise to even stronger radiation in longer wavelength, e.g. the IR, 
as opposed to the UV/optical. We discuss in §2 the pair produc- 
tion, and in §3 the effect of pairs on the prompt flashes from 
reverse shocks. We discuss the current observation of prompt 
optical radiation in §4. Conclusions and discussion on future 
observations are given in §5. 

2. PAIR LOADING IN GRBS 

2.1. GRB emission site 

Consider a GRB central engine producing a relativistic wind 
outflow, which carries energy from the center to some radius 
that GRB emission arises. We assume the isotropic geometry, 
which is valid even for a jetted GRB as long as the Lorentz 
factor r > 1/6 jet, with 6jet the jet opening angle. The width 
of outflow is A < 10^^ cm for a wind lasting a duration of 
< 100 s. The emission energy may come from two energy 
forms, the kinetic energy of the baryon bulk and the magnetic 
energy of the fireball. Therefore the initial wind/fireball en- 
ergy E may be divided into two components (e.g. Zhang & 
Meszaros 2002). The first component is the bulk kinetic en- 
ergy Ek that, as considered in the simplest fireball model, 
comes from the thermal energy of photons, pairs and baryons 
after initial acceleration. The second component is the en- 
ergy Ep of the Poynting flux from the central engine. Define 
a = Ep/Ek, then the fireball is Poynting-flux dominated if 
fj ^ 1, as suggested by recent detection of extreme polariza- 
tion in GRB 021206 (Coburn & Boggs 2003), or kinetic-energy 
dominated if cr < 1. 

The non-thermal spectra of GRBs require that the site of 
GRB emission should be beyond the photospheric radius at 
which the outflow becomes optically thin to scattering by elec- 
trons associated with baryons. The total baryonic electron num- 
ber in the firebaU is Nb = E/{1 + a)VmpC^, with T the 
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Lorentz factor of the fireball. Thus the related optical depth 
Tb = (TrNe/A-KB? <^ 1 leads to a lower limit to the emission 
radius, 



Rh 



3.4xio"4fr3-f (1 



-a) 



-1/2 



cm. 



(1) 



Hereinafter, we use c.g.s. units and the convention Ux = 
U/IQ", except for the Lorentz factors which are scaled as 

— 1/2 

7a; = ^/x. Since Rt (x Tf^ is insensitive to Tt, Rh in the 
above equation do not change much even though the GRB pho- 
tons undergo a few times scattering, T6 ~ a few (not possible to 
be quite large). It is seen that usually i?5 S> A, so the emission 
region of GRBs can be regarded as thin shells. 

GRB light curves at gamma-ray energy are often highly vari- 
able with time scale 6t ~ 1 ms — 1 s. Because of radiation 
beaming an observer can see only the part of the emitting fire- 
ball sphere up to an angle F"^ from the line of sight. The sphere 
curvature leads to a delay in the photon arrival time by R/2r^c. 
This so-called angular spreading time should not be larger than 
the observed variability time scale St, which yields an upper 
limit to R 



Rv 



5.4 X W^^TLM-i cm. 



(2) 



Thus, from the properties of GRB spectra and light curves we 
can constrain the emission site to the regime 



Rb ^ R < Rva 



(3) 



We may take the typical radius as i? = 10^^iii4 cm in the fol- 
lowing. 

2.2. Pair loading 

For an initially thermal-energy-dominated GRB fireball, pair 
production in GRBs begins in the fireball acceleration phase 
when the fireball is still opaque to photons, but the pair number 
density decreases exponentially with increasing radius during 
this phase (Shemi & Piran 1990). When the fireball becomes 
transparent to Compton scattering, the pair number has been 
negligible compared with baryonic electrons. We neglect the 
pair production in this case. 

The intense pair production may occur in the bursting phase 
due to 7— 7 absorption between GRB photons. For the most en- 
ergetic photons at the high-energy end of the spectrum, the opti- 
cal depth of 7 — 7 absorption may exceed unity. Thus, they wiU 
not escape freely from the fireball, and will produce e^pairs. 
The pairs will remain in the fireball, with the same bulk Lorentz 
factor as the fireball (static in the comoving frame). We here 
calculate how many e^pairs generated in this process. 

The observed photon spectra of GRBs can be approximated 
by a broken power-law, with a high-energy portion of the form 
TVe oc for e > £{,, where e = hu/m^c^ is the photon energy 
in units of the electron's rest energy, and 65 ^ 1 is the break en- 
ergy above which the index /3 ~ 2 — 3. The photon number 
above a certain energy e is, therefore, approximated as 



N- 



>e 



P -Ij CbnieC^ 



-(/3-1) 



(4) 



where is the GRB energy emitted in gamma-rays. 

A photon with energy e may annihilate any photons above 
Con = r^/e. Its optical depth due to 7 — 7 absorption is approx- 
imated in an analytical approach by Lithwick & Sari (2001) as 



'77 



(11/180)otA^> 
47ri?2 



(5) 



A photon e with T^j{e) > 1 would be absorbed in the fireball. 
Setting T^j{€cut) = 1 results in a cut-off energy 



ecut — 400£;^_^e(, ^^TIqqR 



14 



(6) 



where the numerical coefficient on the right-hand side corre- 
sponds to typical parameter f3 = 2.2 (Preece et al. 2000). The 
observed GRB spectrum above cut-off energy ecut should be 
attenuated. EGRET on board the Compton-GRO satellite have 
detected several GRBs emitting very high energy gamma-ray 
tails. If the cut-off energy can be detected by future satellites 
such as GLAST, we can constrain the emission radius and the 
initial Lorentz factor, which are the most important physical pa- 
rameters in GRBs. For examples, if no significant attenuation 
above 1 GeV, ecut'meC^ > 1 GeV, this yields a limit 

1477.1/2 ^(/5-2)/2^-(/3-l) 



2.6 X 10^'*S^'-2e6 



.300 



R>RcutilGeY) 

(7) 

The photons above ecut would be absorbed for produc- 
tion, and each absorbed photon corresponds to a pair of e^. 
Integrating the photon spectrum above Ccut we have the total 
number of the resulting pairs 



N± = 1.5 X 10^^i;^_52e, 



2(/3-2)^-2(/3-l) 



■ 300 



i?7 



(8) 



and hence the ratio 



E^;{l+a)R 



14 



(9) 

where we have scaled both E and E^ in units of 10^^ ergs as- 
suming that the efficiency of a GRB is significant. We note that 
the calculation of number in Eq. (8) does not depend on 
detailed models of energy dissipation in emission region, e.g. 
internal shock model or magnetic dissipation model. Our calcu- 
lation is based on observations of GRB properties, i.e. spectral 
and temporal profiles. 

Eq. (9) infers that k± is sensitive to R. The condition for 
strong pair loading, i.e. 2k± > 1, is 



R < Rload = 1.2X10'^ E^,52€t%^^ 



E^,'/\l+ay/\m. 



(10) 

The case 2k± > nip/me in which pair mass dominates baryon 
mass is even achieved if R < Rdom = Rioadi'mp/nie)^^^^. 
From Eqs. (9) and (10) we can see that a wind/fireball with 
typical parameters generally becomes pair rich in the bursting 
phase. These secondary pairs are possible to form an opti- 
cally thick screen again for the GRB emissions (e.g. Guetta, 
Spada & Waxman 2001 ; Kobayashi, Ryde & MacFadyen 2002; 
Meszaros et al. 2002). The optical depth due to Compton scat- 
tering by pairs is t± = aT{2N±)/A'KR^. Substituting Eq. (8) 
into T± = 1 leads to 



R^ = 2.4 X 10^%%ei^-'^^'T^i^-'^^'c^. 



(11) 



A classic GRB takes place in the site optically thin to scattering, 
so the emission radius should be large with R ^ R±. Also the 
permitted regime of Eq. (3) for optically thin GRBs is changed 



3 



to R± <^ R < Rvar- If ^ < the wind/fireball is opti- 
cally thick. However the generated pairs cannot build up an op- 
tical depth substantially exceeding unity, and Compton multi- 
scattering degrades the photons into X-ray energy. This gives 
an interpretation of X-ray flashes (Meszaros et al. 2002)^. 

Fig. 1 shows the above special radii in the F — i? diagrams, 
where the GRB parameter space is separated into different re- 
gions. 

3. PROMPT FLASHES FROM PAIR-RICH REVERSE SHOCKS 

After the bursting phase, the fireball cools into a cold shell, 
which continues to expand at an ultrarelativistic speed and 
with kinetic energy Ea- This afterglow energy should be 
less than the initial wind/fireball energy after the burst, and 
Ea <j E — Ery. When the cold shell runs into the surrounding 
medium, the kinetic energy is dissipated into internal energy, 
which gives rise to prompt emission and long-term afterglow. 
Sari & Piran (1999a, 1999b) and Meszaros & Rees (1999) (see 
also Meszaros , Rees & Papathanassiou 1994; Meszaros & Rees 
1997; Panaitescu & Meszaros 1998) have discussed the prompt 
emission in the context of GRB 990123, assuming no-pair fire- 
ball. We discuss here the prompt emission from a pair-rich fire- 
ball. 

The interaction between the shell and the surrounding 
medium forms two shocks: a reverse shock that propagates into 
the cold shell, increasing its internal energy then giving rise to a 
prompt flash; and a forward shock that propagates into the cold 
medium and produces a long-term multiwavelength afterglow. 
The shocked medium and shell material are in pressure balance 
and separated by a contact discontinuity. There are four regions 
of distinct properties: the unshocked surrounding medium (de- 
noted "1"), the shocked medium ("2"), the shocked shell mate- 
rial ("3") and the unshocked shell ("4"). 

Now the reverse shock may propagate into a pair-rich shell, 
with lepton number density rig, 4 = (1 + 2k±)rn,A and mass 
density p/^ = ni,,4inp[l + 2k±{ine/mp)] (both in rest frame 
of region 4). The shocked material in region 2 and 3 moves 
together, thus they have the same Lorentz factor, 72 = 73. 
The Lorentz factor of the reverse shock relative to region 4 is 
73 4 ~ 74/73. The shock jump conditions yield that the inter- 
nal energy and lepton number density of the shocked shell ma- 
terial are U3 ~ 471 4^40^ and ne,3 — ^"l3Ane,A^ respectively 
(Blandford & McKee 1976). Let a fraction of the internal 
energy goes into leptons, including the baryonic electrons and 
secondary pairs. These shocked leptons are expected to be in a 
power-law distribution: dUe 3 / d7e oc 7"^ for 7e > 7^. with 

The emission reaches a peak when the reverse shock crosses 
the irmer edge of the shell at 



tcr = max[A/c,tdec] 
(Sari & Piran 1999b) with the deceleration time 



, _ o 4 pl/3 -1/3 -8/3 , 
l-dec — 'J-*-C'a.52"-l 74,300 ^- 



(12) 



(13) 



The Lorentz factor of shocked shell at tcr could be estimated as 



73 = 



3^a 



1/8 



= 230 



52 



1/8 



t 



-3/8 
cr.l * 



(14) 



We will focus now on the emission at the peak time tcr- 

The characteristic Lorentz factor (in the rest frame of region 
3) of the accelerated leptons is therefore 



In 



ne,3meC^ 



6^(74/73), 



(15) 



where /j. is the effective mass per lepton in units of electron 
mass and is given by 

linie l + 2k±{mc/mp) 
rup 1 + 2k± 

C 1 for k± < 1; 

~ I (2fc±)-i for 1 <2k±< mp/me; (16) 
[ nic/mp for 2fc-i- > mp/nic. 

Provided that the characteristic Lorentz factor of leptons is not 
less than that of baryons, the value is limited to the regime 



2k±me/mp ^ ^ 1 for 2k± < nip/me, 
~ 1 for 2k± > nip /me- 



(17) 



3.1. The 1 < 2k± < mp/me case 

If k± ^ 1, the random Lorentz factor of leptons becomes 
lower since the lepton number increases significantly. However, 
the pairs do not change the magnetic field in region 3, which 
is assumed to have an energy density times of the internal 
energy density, that is, B = 73c(327rnimp^B)^/^ with m be- 
ing the medium number density. In low-cr case, the magnetic 
field may be shock-induced, thus we assume S,b ~ 10~^, as 
is apparent in the long-term forward shock-produced afterglow. 
But in high-c7 case it may be dominated by the original wind 
magnetic field, we then expect much higher value. Thus, 
for 1 < 2k± < mp/me, the observed peak frequency of syn- 
chrotron photons I'm oc 7^-B is lower by a factor of {2k±)^^. 
Substituting Eqs. (15) and (16) we have 



l^m = 737r, 



eB 
27rm„c 



= 8.2 X 10'"C,-i^ 



10f2 el/2„l/2. 



B 



'±,274,300 



Hz, (18) 



where 73 is cancelled out in the caculation. Note that the emis- 
sion is in quite low frequencies, in contrast with the optical/UV 

for the case without pairs (pm,e/mp ss 1). The cooling fre- 
quency is relevant to those leptons cool at the dynamical time 
t. For t = tcr it is 



1.1 X 10^%'/'^a,52 



E-li\7H:}(^Rz, 



cr,l 



(19) 



most sensitive to ^b- If ~ 10~^, Uc is larger, ~ 10^^ Hz. 

The pairs also significantly increase the radiative efficiency 
by increasing the radiating lepton number. At the crossing time 
tcr, the reverse shock has swept up all the shell leptons, and the 
observed emission reaches a peak. The total lepton number is 
Nc = Nb{l + 2k±) ~ 2N±, and, with help of Eq. (14), the 
peak flux is 



F.^ = 73iVe 



e^B 



AiTdfmeC^ 



32TmimpC^t^j, J \ m 

^Strictly speaking, the way we calculate the pair number is only valid for optically 
no more than an order of magnitude. Furthermore, the calculation of pair nimibers is 
optically thick case here is still valid. 



^1/2 1/4 „ „l/4 , 

' 960 ' — 5-75- Jy, 

dhs^3oo{l + a)t%^ 



(20) 



thin bursts. However, t± is not large, and the spectral break energy is degraded 
insensitive to the break energy provided /3 ~ 2. Therefore, our calculation for 
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where di is the GRB luminosity distance. The flux increases by 
a factor of 2k± as opposed to the case without pair effect. 

The synchrotron spectrum of power-law-distributed leptons 
is simply described by power-law segments: (1) For u < Vm 
we have the synchrotron low-energy tail F,^ oc i^^/^. (2) For 
i^rn < < i^c we have the typical synchrotron slope depending 
on the lepton's distribution, F^, oc (3) For f > Uc 

the spectrum behaves as a fast-cooling slope, F^, oc v~'pI'^ . 
We are particularly interested in the optical prompt emission, 
say, in the R band at vji = 4.5 x lO-"^^ Hz, a frequency well 
above the characteristic frequency and around the cool- 
ing frequency i^c- Because of the independence of Uc on k±, 
no matter whether i/ji is above i^c or not we usually have 



oc 



,(p-l)/2 



OC k '. 



-ip-2) 



Since p k. 2, this means 



that the emission is insensitive to pair loading. For simplicity 
p = 2, we have, for small or vr < v^. 



F, 



0.41- 



1/2 



-^52 -E'y52 74,300 



and for enough large £,b or vr> Vc, 

'fe,-l-£'5274,300 



Jy, (21) 



VR>l'c 



5.8 



cr.l 



Jy, 



(22) 



where k± is cancelled out, and F and 74 can be also cancelled 
out if F ~ 74. 

Synchrotron self-absorption may be important at low fre- 
quencies here. We follow the simple way by Sari & Piran 
(1999b) and Chevalier & Li (2000) to estimate the maximal 
flux emitted by the shocked shell material as a blackbody 



3.2. The 2k± > mp/me case with high a 

A low-fT GRB usually do not produce pairs with mass 
larger than fireball baryons, 2k± > nip/me ^ 1800, because 
that Rdom < Rb (see Fig. 1). However, if cr ^ 1 where the 
initial wind/fireball is dominated by magnetic field, it may ap- 
pear that Rdom > Rb- So, the 2k± > nip/me case may emerge 
below Rdom and above Rb in Fig. 1 for high a. It is usually 
apparent that Rdom < R±^ except for a ^ 10^, therefore the 
case 2k± > mp/me should be associated with optically thick 
pair screens, and, as suggested by Meszaros et al. (2002), ac- 
companied with X-ray flashes^. 

For ~ 1 and ^ w 1 in this case, from Eq. (15) the char- 
acteristic lepton Lorentz factor is simply 7™ « 74/73, which 
is independent of the pair number or k±. The corresponding 
emitted frequency is 



~ 1.1 X 



10^°$^/'n}/%|,,oo Hz, 



(28) 



independent of and k±, where we have scaled 74,100 = 
74/100 for regime below Rdom in Fig- 1- At the same time, 
the Eqs. (19), (20) and (26) above are still vaUd. It should be 

noticed that the peak time tcr is larger for lower Lorentz factor, 

, ^ , _ CQpl/S -1/3 -8/3 

emission, which is dependent on k±, 
k±. It is, for un < Uc, 



We calculate the R-band 



OC F„ u, 



(p-l)/2 



OC 



0.80 



^3/4 
SB, -2' 



1/2, 



^^52^^a!5274,100fc±4 



^L28ri00Cri^cr,2 



Jy, (29) 



F.,66«7r(^) ^kTeff 



(23) 



with the observed size of the fireball R± given roughly by 

RjL « 273C< 
and the effective temperature by 



(24) 



kT, 



eff 



737j.meC^/3, 



(25) 



where 7^ is the Lorentz factor of the electrons that radiate at the 

frequency 1/ and is given by {2'KmeCv / ^fyeBY/'^ . Substituting 
Eqs. (14), (24) and (25) into Eq. (23) leads to 



F„ < F, 



vfib 



150 



5/2 



.-1/4^1/4 



^a,52"'l 



1/2 



-2 .5/4 



7 

-'i,28''cr,l'^y, 

(26) 

where we scale the observed frequency in units of R band fre- 
quency vr. It is obvious that the peak flux in characteristic 
frequency Vm usually suffers strong self-absorption. If the ab- 
sorption frequency Va is defined by Fj, = Fi,^bb, we should have 
t^a> fm- Assuming further < i/c, we can calculate its value. 



2.0 X 10' 



^e,-l^B'^l-£^5274,300 

r3oo(i + CT)i2^_i 



1/3 



Hz. 



(27) 



or for vr > Vc 



6.9 



£^5274,100fc±4 

2griooo'ii cr,2 



Jy, 



(30) 



where we have taken 1 + a a = lOai. Note that k± and tcr 
are not free and are determined by other parameters. 

3.3. Comparision of spectral profiles with fc± <C 1 case 

Regardless of pair effects, the typical frequency of the re- 
verse shock falls in the optical/UV regime, and self-absorption 
can hardly play a role at the optical (e.g. Sari & Piran 1999b). 
So the flux behaves as the synchrotron tail F,^ oc z/^/'' just be- 
low the optical, and Fi, oc at somewhat lower frequencies 
where self-absorption is important. 

For pari-rich reverse shocks, the typical frequency is far 
below the optical and the emission there suffer strong self- 
absorption, so a steep slope 5/2 is apparent in the self- 
absorption regime of lower frequencies. The spectral indices 
from low to high frequencies are in the order of [5/2, ~ 
— 1/2, ~ —1] for i^a < Vc (or [5/2, ^ —1] if Va > /^c)> as shown 
in Fig. 2. We note that the spectral index 5/2 and stronger fluxes 
below the optical are the distinct properties of a pair-rich re- 
verse shock. Future spectroscopy observation of prompt emis- 
sion in longer wavelengths, such as by REM (Zerbi et al. 2001), 
may detect the 5/2 slope in IR spectra, and reveal the pair influ- 
ence. 



This is in the IR band, f „ ~ 4 x 10^^ Hz, for 10"^. As 
shown in Fig. 2, the observed spectrum peaks at Va, and below 
Va it behaves as F^ = F^^bb v^l"^. 

^The generated pair number is prevented from multiplying beyond a value corresponding to r± of a few, due to self-shielding. If t± < 10 is required, the regime 
with R < R±/Vw is not relevant to any observed bursts. So the 2k± > mp/me case, which requires R < Rdom' do not emerge unless Rdom > R±/ViO for 
high a. 
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4. OBSERVATIONS 

ROTSE detected a strong optical flash during GRB 990123, 
which reached a peak of IJy 50 s after the GRB trigger (Ak- 
erlof et al. 1999). From Eqs. (21) and (22) a pair-rich reverse 
shock is easy to emit such a strong flash. Absorption lines in 
the optical afterglow have inferred a redshift of z = 1.6, or a 
luminosity distance of rf; = 4 x 10^^ cm (in a Hq = 71 km 
s^^Mpc^^, il„i = 0.27 and JIa = 0.73 universe). The re- 
leased gamma-ray energy is therefore huge, Ej ~ 10^^ ergs, 
so we assume E ^ Ea ^ 10^^ ergs. Noting that the peak 
time is tcr ~ 50/(1 -f 2) ~ 20 s, the other reasonable pa- 
rameters to achieve a 1-Jy flash are in their typical values: 

-l^|'^i2'^l^%4,300r3oo(l + 0-)"^ ~ 1- 

However, strong optical flashes appear to be rare, since they 
have not been detected from other bursts (Akerlof et al. 2000; 
Kehoe et al. 2001). Especially, GRBs 981121 and 981223 are 
observed to have flashes, if any, fainter than ~ 10 mJy a few 
tens seconds after the GRB triggers. To solve this problem in 
the pair-poor reverse shock model, one needs to assume that 
the characteristic frequency is either well above or below the 
optical, I'm 3> i^R or -c i^r (Kobayashi 2000). However 
Soderberg & Ramirez-Ruiz (2002) have investigated a subset of 
GRBs with afterglow-determined parameters, and shown that 
strong flashes may be expected in all of these bursts. 

For pair-rich reverse shocks, it is usually the case of strong 
flashes in which the optical emission is independent of the pair 
numbers (Eqs. (21) and (22)). Since the gamma-ray fluences 
of GRBs 981121 and 981223 is about one or two orders of 
magnitude lower than GRB 990123, their released gamma-ray 
energy is E^ ~ 10^^ ergs (Kobayashi 2000) if assuming com- 
parable redshift 2; ~ 2 idi^28 ~ 5). To achieve a faint flux, 
< 10 mJy, in Eq. (21) one needs some extreme parameter val- 
ues, e.g. < 10^^. But we expect the fluxes are generally 
large for typical parameters. The rarety of strong flashes in ob- 
servation may need some external reasons to account for, for 
example, dust obscuration of the optical/UV in the star forma- 
tion regions. 

5. CONCLUSIONS AND DISCUSSION 

In this work we calculate the pair number produced in 
GRB fireballs, and study their influence on the reverse-shock 
emission. We have assumed that the intrinsic GRB spectrum, 
before 7—7 absorption, can be extrapolated to very high energy, 
(max ^ ^cut- Our Calculation of pair number does not de- 
pend on any special model, and is only based on observations. 
We show that, for typical parameters, the GRB wind/fireball 
will become pair-rich during the bursting phase, owing to in- 
tense pair production. Later on, the pairs are carried into the 
reverse shock and give rise to the intense prompt flash which 
peaks at lower frequency, i.e. in the IR band, as opposed to the 
optical/UV emission in the case where pair-loading is negligi- 
ble. The emission below the observed peak frequency suffers 
strong self-absorption, leading to the distinct self-absorption 
spectral index 5/2. 



Since k± oc F^'-^^^^'^-'i?^^, pair loading is unimportant for 
Lorentz factors large enough (also see Fig. 1). However, the 
Lorentz factors of GRBs have been found not to be quite large 
by some previous works. For examples, in the internal shock 
model of GRBs, the burst emission should come from internal 
shocks before significant deceleration due to the swept-up ex- 
ternal matter, which results in F ^ 10^ (Lazzaiti, Ghisellini & 
Celotti 1999); The observed pulse width evolution in the burst- 
ing phase seems to rule out the GRB models with F » 10^ 
(Ramirez-Ruiz & Fenimore 2000); For GRBs to be efficient, 
the cooling timescale of leptons should be smaller than the 
GRB duration, leading to an upper hmit of ^ 1200 (Derishev, 
Kocharovsky & Kocharovsky 2001). We, therefore, expect that 
most GRBs are pair-rich with 2k± > 1. 

The pair-rich reverse shocks emit mainly in the IR band. If 
GRBs are associated with star forming regions, as imphed by 
the observations (e.g.. Bloom et al. 1999; Fruchter et al. 1999; 
Bloom, Kulkarni & Djorgovski 2001) and by the favored pro- 
genitor models that GRBs come from explosions of massive 
stars (Woosley 1993; Paczyriski 1998; Vietri & Stefla 1999; 
MacFadyen & Woosley 1999), the optical may suffer dust ob- 
scuration, whereas the IR would not (e.g., Draine & Hao 2002; 
Reichart & Yost 2002; Ramirez-Ruiz et al. 2002; Fruchter et 
al. 2001; Waxman & Draine 2000). Strong pair-loading also 
increase significantly the fiash brightness. Thus, we expect that 
there are many strong IR flashes associated with GRBs (and 
XRFs). If so, the rarity of optical flashes in GRBs will in- 
fer dust-enshrouded environment of GRBs. The detection of 
prompt IR flashes requires rapid response to the GRB trigger in 
such band, and the up-coming near-IR telescope, REM (Zerbi 
et al. 2001), is capable of catching them. 

The observation in high energy range, e.g., ^ 100 MeV, 
plays an important role in GRB research, because they pro- 
vide a way to determine the initial Lorentz factor F (e.g. Bar- 
ing 2000; Lithwick & Sari 2001, and references therein). If a 
GRB exhibits attenuation above some ecut< we can use Eq.(6) 
together with Eq. (2) to constrain F. If Ccut ~ 1 GeV, then 
we have, for typical parameters, F ^ 230 where the equality is 
for R Rvar (scc the crossing point between curves Rcut and 
Rvar in Fig- !)• The high-resolution spectra of GRBs may be 
provided by the next generation high-energy gamma-ray obser- 
vatory, e.g., GLAST. 

Eq.(8) indicates that ~ lO^^'E'^ 52 pairs are produced in a 
GRB. This number should be reduced by a factor of ffjet/'^ if 
GRBs are beamed. It is recently suggested that the annihilation 
of the relic positrons provides identification beacon for at least 
one GRB remnant in the Milky way (Furlanetto & Loeb 2002). 
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Fig. 1 . — Emission site R versus initial Lorentz factor F: <t = (upper frame) and <t = 10 (lower frame). Neitlier tlie regime above Rvar nor tliat below Ri, is 
relevant to observed GRBs. Below Rioad is the regime where pair loading is important, 2k± > 1. Case 2k± > rrip/nie is even available below Rdom- A classic 
GRB, with non-thermal sub-MeV radiation and rapid variabilities, should be optically thin to scattering of photons by baryonic electrons (R > Rf,) and secondary 
electron-positron pairs (R > R±). The X-ray flashes are possible to appear below -R±. The GRBs shows no attenuation above 1 GeV correspond to regime beyond 
the curve marked -Rcut(l GeV). £52 = £'7,52 = ^6 = St-i = 1 and /B = 2.2 are assumed. The scaling laws of special radii with the parameters are also marked, 
except for tf, to which the relation is usually weak. 
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Fig. 2. — Examples of prompt-emission spectra from pair-rich reverse shocks: (^s, ct) = (0.01, 0) {solid line), (0.2, 5) (dashed line) and (1, 10) {dotted line), 
respectively. Other typical parameter values are assumed (see the text), especially that fe±,2 = tcr,i = 1- 



